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It is indicated that there is a possibility of spin polarization in the quark matter at high 
density due to the tensor-type four-point interaction between quarks. The realized phase 
is determined by the comparison of pressure of quark spin polarized phase with that of the 
two-flavor color superconducting phase which may appear in high density quark matter. It 
is shown that the quark spin polarized phase is realized at rather high density, while the 
two-flavor color superconducting phase is realized at a lower density region. 
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§1. Introduction 



Under extreme conditions such as high temperature and/or high baryon den- 
sity, it is interesting to study the behavior of the quark and gluonic matter and/or 
hadronic matter governed by quantum chromodynamics (QCD) in the context of the 
physics of the relativistic heavy ion collisions and of the interior of compact stars, 
psj \ At low temperature and high density region 1 ) such as the interior of a neutron star, 

it is widely investigated and pointed out that various phases may be realized, for 
example, various meson condensed phases in hadronic matter, 2 ) the two-flavor color 
superconducting (2SC) and the color-flavor locked (CFL) phases, 3 )~ 5 ) the quarkyonic 
phase 6 ) and the quark ferromagnetic phase 7 ) in quark matter, and so forth. Espe- 
cially, we may conjecture that the quark ferromagnetic phase exists at high density 
quark matter, since the existence of magnetar was reported. 8 ) In a preceding work, 7 ) 
the pseudovector-type interaction between quarks was considered and it was pointed 
out that the quark spin polarization was realized. However, it was shown that, the 
quark spin alignment disappears, if the quark mass is zero, for example, in the chiral 
symmetric phase. 9 ) 

Recently, the present authors have indicated that the quark spin polarization 
occurs at high density quark matter even in the chiral symmetric phase due to the 
tensor-type four-point interaction between quarks, which leads to the quark ferro- 
magnetization. 10 )' 11 ) In that paper, it has been shown that a second-order phase 
transition occurs from normal quark phase to quark spin polarized phase. However, 
in quark matter at high baryon density, it is believed that the two-flavor color su- 
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perconducting phase appears. Thus, it is interesting to investigate the stability of 
the quark spin polarized phase against the 2SC phase. 

In this paper, we investigate the color-symmetric superconducting phase 12 -* and 
the quark spin polarized phase and clarify which phase is stable by using the Nambu- 
Jona-Lasinio (NJL) model 13 ) including quark-pair interaction 14 ) and the tensor-type 
four-point interaction. 10 )> 1:L ) The quark-pair interaction is derived by the Fierz trans- 
formation from the original NJL model Lagrangian. We consider the quark-pairing 
gap A and the quark spin polarization F separately such that A ^ or F 7^ leads 
to the 2SC phase or quark spin polarized phase, respectively. In this paper, the sys- 
tem is treated by the mean field approximation and the BCS state is introduced. 15 ) 

This paper is organized as follows: In the next section, the basic Hamiltonian 
is introduced and expressions of the Hamiltonian under some basis sets are given. 
In Appendix A, one of the expressions of Hamiltonian is given in which the good 
helicity states are used. In §3, the BCS state is defined and the thermodynamic 
potential is derived. In §4, numerical results and discussions are given and the last 
section is devoted to a summary and concluding remarks. 

§2. NJL model with tensor interaction 

2.1. Basic Hamiltonian 

Let us start with the following Lagrangian density: 

C = Cq + Cs + Ct + C c , 

C s = G s {bH>? + (V^srVO 2 ) , 

Q 

C T = - — ((ip7 tJ 'l u T'ip)(ip7^vTip) + (V^757 M 7*V) (^757^7^)) , 

£c= T £ {{^T2\A^ C )^ C i^T 2 \ A ^) + {^T 2 \ A ^ C )^ C T 2 \ A ^ .(2-1) 
A=2,5,7 

Here, i\f = C$ T with C = ij 2 j° being the charge conjugation operator. Also, t 2 is 
the second component of the Pauli matrices representing the isospin sn(2)-generator 
and X A are the antisymmetric Gell-Mann matrices representing the color su(3) c - 
generator. As is well known, the Lagrangian density Co + Cs corresponds to the 
original NJL model. We add C c which can be derived by the Fierz transformation 
from Cs and represents the quark-pair interaction. Also, we introduce Ct which 
represents the tensor-type interaction between quarks. 

In this paper, we concentrate on quark matter at high baryon density where the 
chiral symmetry is restored in the density region considered here. Thus, the chiral 
condensate, (ipip), is assumed to be equal to zero in this high density region. Under 
the mean field approximation, the above Lagrangian density is recast into 

C MF = C + C T 4F + Cf F , 
C MF = -F^E 3 t 3 ^) - g , 
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F = -G$£ 3 t^) , U 3 = -i 7 V = ( ^ 



3A 2 



2 ^ v r ,o * ~ . / 2G , 

A=2,5,7 C 

Z\ A = ^ = -G c (^7 5 r 2 A A V') , Z\ = Z\ 2 = ^5 = ^7 , (2-2) 



where /i.e. represents the Hermitian conjugate term of the preceding one. Here, we 
used Dirac representation for the Dirac gamma matrices and a 3 represents the third 
component of the 2x2 Pauli spin matrices. The symbol (• • • ) represents the expec- 
tation value with respect to a vacuum state. The expectation value F corresponds 
to the order parameter of the spin alignment which leads to quark ferromagneti- 
zation. The expectation value A corresponds to the quark-pair condensate which 
means the existence of the color superconducting phase if A ^ 0. Here, in order 
to insure color symmetry we assume that all quark-pair condensates have the same 
expectation values, A 2 = A$ = A-j. 

The mean field Hamiltonian density with quark chemical potential \i is easily 
obtained as 

Umf -fiAT = JC + H^ F + Uf F , 
= V>(-7 • V - ^To)^ , 

J^MF _ _qMF H MF — —jd MF (2-3) 



with Af = ip^ip. In the basis of good helicity states, the relevant combination of the 
mean field Hamiltonian Hmf = / (PxHmf and the quark number N = f d 3 xJ\f is 
given by 



H M F -llN = £ [iP~ ^) c Ura C PVr a ~{P + ^CpryraCprrr 
pr/ra 

+F V<A [ vjj+g ( c \ c +c ~t ~ c \ 

pr]ra L 

P3 ( f , ~f \ 

V \^pr S T(x c pr\ Ta ' c prira c PV ra J 

+ "2~ ( C PV raC -priT'a' ^PVTofi-pryr'a' 

prjaa' 'a'' 'tt' 

F 2 3A 2 

+V -2G + V -2G C ' (2 - 4 » 

where V represents the volume in the box normalization.*) Here, c PriTa and c PriTa 
represent the quark and antiquark creation operators with momentum p, helicity 



*-* We hope that no confusion may occur although notations V C3 and V(p) appear later which 
represent the quark-pair interaction in (2-5) and a matrix in (2-9) diagonalizing a certain part of 
the Hamiltonian matrix, respectively. 
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rj = ±, isospin index r = ± and color a. Further, </> T = 1 for r = 1 (up quark) and 
cf) T = — 1 for r = — 1 (down quark). Also, e TT / and e aa ' a " represent the complete 
antisymmetric tensor for the isospin and color indices. We define p = \/p 2 + p\ + V% 
that is the magnitude of momentum. 

For later convenience, we rewrite the above Hamiltonian as 



Hmf — l^N = Hff F 



F 2 3A 2 
,N + V CS + V.- + V. 



( c. 



ttSP 
n MF 



UN 



E 



.t 

r 



\ 



( 



(k - n ■ l)5 TT >5 a 



/ p \ 



Cp+r'a' ^ 
Cp—r'a' 
Cp+r'a' 
\ Gp—r'a' j 



p 
-p 

V o o o -p J 

A 
2 



+ 



Fr 3 



/ _0 VPl+P2 -P3 

VpT+pI o o 

-P3 _0 

o p 3 Vp'i+pi 





P3 

Vp'i+pI 





— 2_^i \ C pr)Ta C -pr]T'a' + C pr)Ta C -p-q T ' a 1 + c -pv r ' a ' c pr\Ta + C-pryr' a' c prrra) 



pr/aa'a"TT' 



As for V cs , we can easily sum up with respect to r and r' which leads to 



(2-5) 



t //(c 



t 



a"\^pr)+a^-pri-a 



i + C 



/ ~i~ C—pri—a'Cpri+a ~\~ C—pri—a'Cpri+a) ■ 

(2-6) 



2.2. Expressions of the Hamiltonian under other basis sets 

If the quark-pair condensate A is equal to zero, the operator given by equation 
(2-5) becomes Hf[ F — fiN + VF 2 /(2G). Since the quark number N is already di- 
agonal, let us diagonalize the Hamiltonian matrix k first. First, we set up T3 = 1 
because the contribution of T3 = — 1 is the same as the contribution of T3 = 1 as was 
seen in (2-6), which results in a factor 2. Then, the Hamiltonian matrix k can be 
simply expressed as 



k = 



( q e -g \ 

e q g 

—g — q e 

V g e -q j 



(2-7) 



where q = p, e = F^Jp\ -\-p\jp and g = Fp^/p. The eigenvalues of k are easily 
obtained as 

±4 ±} = ±Vfl 2 + (e±<7) 2 = ±^pl+(F±Jp 2 +pij . (2-8) 
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The following matrix is introduced in order to diagonalize the above Hamiltonian 
matrix: 



V(p) 



^/4 +) + e +<? 


\Je p ) -e+q 


je^-e-q 


\/4 ] +e~q \ 




2^ 


2y/4^ 


2^ 


\Je p +) +e+q 




\Je p +) -e-q 




9 


9 




2^ 
9 


2 v / 4 +) (4 +) + e +9) 

9 


2^e p -\e p - ) -e+q) 
9 


2 V / 4 +) (4 +) - e -?) 

9 


2^e p -\e p - ) +e-q) 
9 


2^4 +) (4 +) + e + (? ) 


2^ e ( p ] (e ( p } -e+q) 


2 V / 4 +) (4 +) - e -9) 


2 s Je p -\e ( p - ) +e-q) ) 



V 

Then, we can diagonalize the Hamiltonian matrix k as V*kV. Namely, 

\ 



(2-9) 



ttSP 
n MF 



£ 

pra 



\ p—TCt 



diag k 



bp— ret 
bp+ra 



(2-10) 



where we define 



diag k = 

bp— ret 
bp-\-ra 

\ bp— TOl 



.(+) 
-p 







-(-) 

-p 










— e 



(+) 









V\p) 



J 



^p—ra 
Cp+ra 



(2-11) 



Secondly, let us consider the total Hamiltonian matrix with quark-pair interac- 
tion term in the above derived basis. In this basis, V cs in (2-6) is written as 



V CS = A 



u p++ct 
u p — ha 
U p++a 

v ) 



vHp)v(-p) 



( h -p + - a > \ 

—p a' 



(2-12) 



Here, if we replace p into — p, then q, e, g and in (2-7) are changed or unchanged 
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to q, e, —g and Thus, we obtain 



/ q+e 
1 ,(+) 



V\p)V{-p) 



,(+) 



V 





g-e 




_JsL n \ 
,(+) u » 





g+e 




,(-) 



"4^ / 



(2-13) 



Finally, let us diagonalize the above matrix V\p)V{— p). We introduce the 
following matrix: 



/ 



W 



4 +) + e +g 



2s. 



\g\ 



T+T 







4 -e+g 



^24 +) (4 +) +e+ g ) 

o 



2e 





(=) 



2 £ _ 



\g\ 



T+T 







^24 +) (4 +) -e-g) 





2s r 



Iffl 



FT 



1/24 '(4 ) + e ~9) / 

(2-14) 



Then, by using the matrix the matrix W(p)F(— p) which appears in the quark- 
pair interaction part can be diagonalized as 



wW\p)v(-p)w 



By introducing new fermion operators (ttp^ra? ^prjrai dpr/raj QprjTct) by 



/ 1 











\ 





1 
















-1 







V 








-1 


/ 



(2-15) 



U p—Ta 

a p+ 

TOt 
\ — TOL 



bp— ret 
bp-\-ra 

\ bp—TU ) 



(2-16) 



the mean field Hamiltonian in which both the quark spin polarization and quark-pair 
condensate are simultaneously considered can be expressed as 



pr]TOL 

A 



+- 



a) a 1 " - - f 

p-qra' —p-qr'a" 



p-qra^priTa 
t 



a' ,a' „ + h.c. 

prjTCt' ~pr)TCt" 1 



e aa ' a "e TT Kp T 



F 2 3A 2 
+V -2G +V -2G C 



(2-17) 
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Here, the summation with respect to isospin indices r and r' is explicitly rewritten. 
The above mean field Hamiltonian is the starting point for discussing the quark spin 
polarized phase and color superconducting phase. 

§3. BCS state and thermodynamic potential 

We introduce the BCS state following Ref.15): 



S= E ¥ E 



2 a pr]Ta a - prl T'a' e aa'a"^TT'YT 



+ ^ ] - ^ ] 0'pr]raO'-prjT'a' e aa'a" e TT' ( t > T j 

where |0) is the vacuum state with respect to 

i^_ pr; are satisfied. Here, the contribution of quark-pairing with negative energy is 
not considered in quark matter. Then, the state \W) is a vacuum state with respect 
to new operators d pr)Ta : 

J a priTa - K pf] (a} , a - a) , ) for e$ > /x . . 

a pV ra ~ \ ] ,r> ( - \ f fa) <r 

where (a, /?, 7) is a cyclic permutation for color indices and r' = — r. For the above 
operators, the state |#) satisfies the following relation: 

d^W) = . (3-3) 

The state is identical with the BCS state. 

The expectation values for the BCS state are summarized as follows: 

X PV = (^l a pr,™ a pr,™'|^> , ( for OL^o!) 

Np,) = { , ^ r \cip 11TO ,ap 1 j Ta \W r ) , 

D pv = {^\a^ p ^ Tal a mTa \^) , (for a / a') 

m , (3-4) 

where we take a = 1 and a' = 2 for the following calculations because the color 
symmetry is retained. We can easily calculate the above expectation values by using 
the relation (3-3) with (3-2). 

(v) 

For £ p > fi, the following relations are obtained: 

X pr) = -Kp, q + Kp V (N PV - X m ) , 
Np,j = 2K pfl — 2Kp fl (N prj — X pr) ) , 

Djyq = K pr i — 3K prj D pr i + K prj P prj , 

P p -q = —^Kp V P pri . (3-5) 
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(v) 

In the same way, for e p < ji, the following relations are also obtained: 



X 



PV 



pq 



l + 2K 2 pr) (X PV -N PV ) , 



D — K — 3K 2 D — K 2 P 

L -'pr\ — ly pq oli -pij ly PV 1 ^prj Jr PV ! 

P — 2K 2 P 



(3-6) 



As a result, the expectation values are calculated as 
For e p v) > /j, 



K 2 

Y — _ P T > 

m ~ 1 + 3K 2 



PV 



_ 2K 2 V 
' V '" 1 + 3K 2 



pn 



1 + 3K 2 



P PV = , (3-7) 



PV 



For e p v) < n 



Xprj — 



K 2 



1 + 3K 2 V 



D 



PV 



PV 



1 + 3K 2 m 



PV 



PV 



1 + 3K 2 pr] 



. 



IK 2 
1 + 3K 2 V 



(3-8) 



Thus, by summing up with respect to color and isospin indices, we obtain the ex- 
pectation value of the mean field Hamiltonian for the BCS state as 



(9\H MF -^N\9) = 3 Yl 

pr?r(4'' ) >M) 
PVrie^ <fi) 



2K 2 
A*) m 



1 + 3K 2 m 



K 

+ 2A 



1 + 3K 2 prj 



(4" } 



2K 2 
l + 3K 2 ^ 



+ 2A- 



K, 



PV 



l + 3K 2 m 



F 2 3A 2 

+v -2g + v -2g: 



(3-9) 



In order to determine the BCS state, namely, to determine the variational variables 
K^pfj ehicL K^pyj j we introduce new variational variables 9 prj and 9 PV instead of K PV and 



PV 



sin 6. 



\[3K, 



PV 



sin i 



PV 



PV 



PV 



cos 9 



V 1 + 3K pv 



COS I 



PV 



PV 



V 1 + 3 ^ 



(3-10) 



Thus, the Hamiltonian and color superconducting gap A can be expressed in terms 



of 9 pri and 9 prj as 



(&\H M f - nN\&) = 2 ^(ej,^ - A*) sin2 ^ + 2>/3^ sin 6> pr , cos pr , 



pv(^>h) 
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+2 ^2 [(e^ - fi) (3 - 2 sin 2 6 PV ) + 2y/SA sin 6> pr/ cos 6 



F 2 3Z\ 2 
+y -2G +y -2G: 



(3-H) 



Z\ = A 2 = A 5 = A 7 = -G c • i £ (^la^a 



-prj—T'y 



pri/3-yT 



-2G C . 



4G C 



y 4^ 1 + 3^, k 4^ 1 



K. 



PV 



- y 

v ^ 



sin t/ PV cos 



PV 



- y 

V ^ 

PV{£ { p ] <^) 



sin t/ pr) cos 



(3-12) 



Here, in the first line in (3-11) and the third line in (3-12), the isospin indices are 
summed up and an extra factor 2 appears with respect to the corresponding expres- 
sion in 12). Also, in the second line in (3-12), color indices, j3 and 7, are summed up 
and an extra factor 2 also appears. 

Next, we impose the minimization condition for (\P\Hmf ~ nN\\P) with respect 
to A, 9 pv and 9 prj . First, we impose the minimization condition with respect to A: 



d_ 
OA 



(<P\H M F - fiN\&) = . 



(3-13) 



This minimization condition leads to the gap equation in (3-12) exactly. Secondly, 
we minimize (\P\Hmf — ^N\^) with respect to 6 prj and 6 prj : 



-{V\H MF - fiN\&) = , 







(V\H MF - fiN\$) = 



(3-14) 



which lead to the following equations: 

(Ep^ — fi) sin 26 PV + y/3A cos 26 PV = , 
-(4^ - aO sin 26 PV + V3A cos 29 PV = . 



(3-15) 



Thus, we obtain tan 26 PV 
gives 



y/3A/(e$P - fj) and tan2f9 pr/ = y/%A/{e^ } - n) which 



fa) 



sin 1 



sini 



pr; 



PV 




^ (e p v) - fi) 2 + 3A\ 



(3T6) 
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Since the variational parameters 6 pri and 9 PV are determined completely, namely, K pv 
and K pr) are determined, the BCS state is obtained. 

Thus, we can derive the thermodynamic potential <P(A, F, jj) at zero temperature 
from (3-11) with (3-16) as 

${A,F,n) = ±-(<p\H MF ~ ^N\$) 



= 2-1 £ [2(4 T,) - M) - y/ieP ~ ») 2 + 3^ 2 

1 A 

«4 E 



V(4" j -m) 2 + 3^ 2 



fa) 



F 2 3Z\ 2 
+ 2G + 2G~ C 



(3-17) 



where e 



(±) 



/pg + ^-F ± v / Pi~+~pI) • Here, we explicitly introduce and write a 

three momentum cutoff parameter A The gap equation in (3T2), that is d<P(A, F, /j,)/dA 
0, is obtained as 



1 A 

2-1 y 

V ^ 



1 

G~c 



§4. 



pr?=± -//) 2 + 34 2 

Numerical results and discussions 







(3-18) 



In this section, we investigate the possible phases with F = and A ^ or with 
F / and A = at high density, namely, two-flavor color superconducting phase 
or spin polarized phase, respectively. First, in the quark spin polarized phase with 
A = 0, the thermodynamic potential reduces to 

#(4 = 0,^) = 6-1 £ (4 +) -p) + 6-l £ ( e H + (4-1) 



p (4 +) <m) 



p (4 '<m) 



where it should be noted that \J (e^ — /u) 2 = £p 7; — for > p or — (£p 7; — //) 

for < fi, respectively, when (4-1) is derived from (3T7). Of course, the sum over 
momentum is replaced by the momentum-integration as 



_(*?) 



(v) 



d 3 p 
(2vr) 3 



(4-2) 



Then, the thermodynamic potential derived here is identical with one which has been 
previously obtained by the present authors, that is, Eq.(3.1) in Ref. 10).*) Thus, the 



In Eq.(3.1) in Ref.10), F-integration has to be carried out. 
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thermodynamic potential with A = has already been given in 10) in the analytical 
form as 



F^_l_ 
2G~^ 



(3F 2 n + 2/j, 3 ) +F/x 3 arctan 



0(A = Q,F,n) = { 



"T ln 



fi+Vv 2 - F 2 



F 2 1 ,^ 

ri 3 F 

< 2G 2vr^ 



for F < n 

for F > n 

(4-3) 



If the minimum of thermodynamic potential exists in the range of F < fi, the gap 
equation for F is derived from d<P(A = 0, F, jj)/dF = which leads to 



F=C 

IT 2 



2Ffiy / n 2 -F 2 + // 3 arctan 



3l ^ + v^^ 



V/" 2 - F 2 



F 6 In 



(4-4) 



In the other case, namely F = case, the two-flavor color superconducting 
phase may be realized with A ^ 0. When F = 0, the quasiparticle energy = p 
is obtained for rj = ±. The thermodynamic potential (3-17) is then evaluated as 



${A, F = 0, fi) = 4 • i ^ (2(p - M ) - V(p-/^) 2 + 3^ 2 ' 

+4'-^((p-^)- V(P " + 3^ 2 ) + 



3Z^ 
2G r 



(4-5) 



After replacement of the sum over momentum into momentum-integration in (4-2), 
the analytical form of the thermodynamic potential is obtained with a three-momentum 
cutoff parameter A: 



$(A,F = 0,») 



3^__^_ A 4 2fiA 3 
2G~ C ~ 6^ + 2k 2 ~ 3vr 2 



1 



127T 2 



(2^ 3 - 39A 2 ii)(y / f i 2 + 3A 2 - {^ - A) 2 + ZA 2 ) 
+ (6/l 3 + 9A 2 A - 2A 2 n - 2Afi 2 )^(fi- A) 2 + 3A 2 
+ 3(12AV-9A*)ln A -» + ^- A y 2 + 3A2 



-H + yV + 3A 2 



(4-6) 



The gap equation becomes 



1 



7T^ 



-3/zyV + 3Z^ 2 + (3/i + /l)V(/"-^) 2 + 3^ 2 



+ (2,u 2 - 3zA 2 )ln 



A - ii + y/(fi - A) 2 + 3A 2 



1 



(4-7) 
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Table I. Parameter set 



A 1 GeV 


G 1 GeV~ 2 


G c / GeV" 2 


0.631 


20.0 


6.6 



In both the cases, the quark number density p is calculated from the thermodynam- 
ical relation as 

In our model Hamiltonian, the model parameters are G, G c and A. In the original 
NJL model, the coupling strength Gs appears where Gs = 5.5 GeV" 2 is adopted 16 ) 
although this parameter does not appear explicitly in the model considered here. 
Then, following Ref.14), the strength of the quark-pair interaction, G c , is taken*) as 
(G c /2)/Gs = 0.6. Thus, we adopt G c = 6.6 GeV -2 . As for the three-momentum 
cutoff A, a standard value 16 ) is adopted here, namely, A = 0.631 GeV. If the tail 
stretching beyond the Fermi momentum of occupation number should be fully taken 
into account in the BCS theory, a larger value of three-momentum cutoff should be 
adopted such as A = 0.8 GeV. 15 ) However, it will be later seen that the tail of the 
occupation number is already sufficiently taken into account in the case A = 0.631 
GeV. As for G in the strength of the tensor-type interaction between quarks, we 
put G = 20 GeV -2 , which was used in our previous paper. 10 ) As was discussed in 
Ref.10), if the effect of the vacuum polarization is taken into account, the coupling 
constant G should be replaced to the renormalized coupling G r in which 1/G r = 
l/G — A 2 /tt 2 . Then, G r = 20 GeV -2 corresponds to the bare coupling G = 11.1 
GeV -2 for A = 0.631 GeV or G = 8.7 GeV -2 for A = 0.8 GeV. Thus, the strengths 
of the quark-pairing and tensor-type interactions are comparable. The parameter 
set used here is summarized in Table 1. 

The thermodynamic potential is numerically estimated in two phases, namely, 
color superconducting phase with F = and A ^ and quark spin polarized 
phase with F / and A = 0, separately. We compare the pressure in the color 
superconducting phase with that in the quark spin polarized phase. The pressure p 
is given by 

p = -$(A,F, f i) . (4-9) 

In Fig. 1 (a) , the pressures for normal (thin curve) , two- flavor color superconducting 
(dash-dotted curve) and quark spin polarized (solid curve) phases are shown in the 
case of A = 0.631 GeV as functions with respect to quark chemical potential fi. Up 
to n = 0.442 GeV, the two-flavor color superconducting (2SC) phase is realized. 
However, above \i = 0.442 GeV, the quark spin polarized phase is realized. In Fig. 
1 (b), details are depicted around /j, ~ 0.442 GeV. The value of chemical potential 
in the phase transition point is 0.442 GeV, which leads to the 2SC phase up to 
H = 0.442 GeV. Thus, it is enough to include the effects of the tail of occupation 
number sufficiently with three-momentum cutoff A = 0.631 GeV. 



In Ref.14), Gc/Gs = 0.6 is adopted where Gc in Ref.14) corresponds to G c /2 here. 



13 




0.442 WlGeV] 0.442 UlGeV] 

(a) (b) 

Fig. 1. (a) The pressures for normal (thin curve), two-flavor color superconducting (2SC) (dash- 
dotted curve) and quark spin polarized (SP) (solid curve) phases are shown as functions with 
respect to quark chemical potential n in the case of A = 0.631 GeV . (b) The details are depicted 
around y, w 0.442 GeV. 



Since we are interested in the high density region where chiral symmetry is 
restored, we focus on the transition point from color superconducting phase to the 
quark spin polarized phase. In Fig. 2, the realized phase is shown by the solid 
curve. The vertical axis represents the pressure and the horizontal axis represents 
the baryon number density divided by the normal nuclear density po = 0.17 fm -3 
instead of the quark chemical potential by using the relation (4-8). Up to ps = 
4.73poj the color superconducting phase is realized. From ps = 4.73po to 5.84po, the 
color superconducting phase and the quark ferromagnetic phase may coexist. Above 
Pb = 5.84po, the quark spin polarized phase is realized. Thus, the phase transition 
may be the first order one. 
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Fig. 2. The realized phase is shown by the solid curve in the case of A = 0.631 GeV. The vertical 
axis represents the pressure and the horizontal axis represents the baryon number density divided 
by the normal nuclear density po = 0.17 fm -3 . 
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§5. Summary and concluding remarks 

In this paper, it was shown within the mean field approximation that the quark 
spin polarized phase is realized succeeding to the two-flavor color superconducting 
phase as the baryon density increases. We calculated the pressure of the two phases 
separately, namely, in the 2SC (A ^ and F = 0) and the quark spin polarized 
(A = and F^O) phases. As a result, at a certain density region, the 2SC phase is 
realized. However, at rather high density region, the quark spin polarized phase is 
realized. If the large F expansion is carried out in Eq.(3T8), then the gap equation 
for A is expressed approximately as 

\ f d 3 p ( 2»^p{Tp z 2 -3A 2 
^ [J (2vr)3 ^ + 2F2 

For large F, the above equation has only the trivial solution, A = 0. Thus, it is 
expected actually that the quark spin polarized phase with A = and F ^ is 
realized in the high baryon density region with large F. 

For three-momentum cutoff A = 0.631 GeV, the tail of the occupation number 
is sufficiently taken into account since the value of chemical potential on the phase 
transition point from 2SC phase to the spin polarized phase is 0.442 GeV. In general 
in the system with finite chemical potential, a chemical-potential dependent cutoff 
may be used. 17 ) However, for the sake of comparison, if a rather large fixed value of 
three-momentum cutoff, A = 0.8 GeV, is adopted where the tail of occupation num- 
ber is more fully taken into account, the critical chemical potential from 2SC phase 
to the spin polarized phase is obtained as p, = 0.491 GeV. Above p = 0.491 GeV, 
the quark spin polarized phase is realized. With this larger fixed three-momentum 
cutoff, the following results are obtained: Up to ps = 7.42po, the color supercon- 
ducting phase is realized. From ps = 7A2po to 10.86po> the color superconducting 
phase and the quark spin polarized phase may coexist. Above pb = 10.86/90) the 
quark spin polarized phase is realized. The use of large three-momentum cutoff leads 
to the shift of the phase transition density, that is, the phase transition density is 
higher with larger three-momentum cutoff. 

One of the next task may be to consider the quark-pairing gap A and the spin 
polarization F simultaneously. Under this treatment, the behavior of phase transi- 
tion from 2SC phase with A ^ and F = to the quark spin polarized phase with 
A = and F ^ may be clear because the coexistence phase with A ^ and F ^ 
may appear. This is a future problem. 

Also, it is important to investigate the interplay between color-flavor locked 
(CFL) phase and the spin polarized phase at higher density region in quark matter. 
Further, it may be interesting to investigate the interplay between two quark spin 
polarized phases, namely, the quark spin polarized phase originated from the tensor- 
type four-point interaction and that originated from the pseudovector-type four-point 
interaction between quarks. These topics are left for future problems. 
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Appendix A 

Mean field Hamiltonian in the basis of good helicity states 

In the Dirac representation for the Dirac gamma matrices, the Hamiltonian 
matrix of the spin polarization part, H^f F , is written as 

hf[ F = p-a + Fr 3 f3U 3 

( Ft 3 
-Ft 3 

P3 Pi ~ iP2 

\ Pi + iP2 ~P3 

where a 1 = 7°7* and j3 = 7 . For good helicity states, this Hamiltonian matrix is 
diagonalized with F = 0. For simplicity, we rotate around ^3 axis and we set P2 = 



P3 
Pi + iP2 
-Ft 3 




Pi ~ iP2 
P3 



Ft 3 



(Al) 
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without loss of generality. In this case, we derive k in §2 as follows: 



U = 



Vp + P3 



Vp -ps 
-^Vp~+p~3 

VP-P3 



-Vp + ps 
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P3 


bil 





/ 



(A-2) 



Finally, in the original basis rotated around p3-axis, \p\\ is replaced to \Jp\ + p\ and 
we get the result given in Eq. (2-5) in §2. 



